Approximation numbers of Sobolev embeddings—Sharp constants and tractability

We investigate optimal linear approximations (approximation numbers) in the context of periodic Sobolev spaces H^s(T^d)$H^s\left({\mathbb{T}}^d\right)$ of fractional smoothness s>0$s>0$ for various equivalent norms including the classical one. The error is always measured in _L2(T^d)$L_2\left({\mathbb{T}}^d\right)$. Particular emphasis is given to the dependence of all constants on the dimension d$d$. We capture the exact decay rate in n$n$ and the exact decay order of the constants with respect to d$d$, which is in fact polynomial. As a consequence we observe that none of our considered approximation problems suffers from the curse of dimensionality. Surprisingly, the square integrability of all weak derivatives up to order three (classical Sobolev norm) guarantees weak tractability of the associated multivariate approximation problem.     

Strong tractability of multivariate integration using quasi-Monte Carlo algorithms

We study quasi-Monte Carlo algorithms based on low discrepancy sequences for multivariate integration. We consider the problem of how the minimal number of function evaluations needed to reduce the worst-case error from its initial error by a factor of depends on and the dimension . Strong tractability means that it does not depend on and is bounded by a polynomial in . The least possible value of the power of is called the -exponent of strong tractability. Sloan and Wozniakowski established a necessary and sufficient condition of strong tractability in weighted Sobolev spaces, and showed that the -exponent of strong tractability is between 1 and 2. However, their proof is not constructive.

In this paper we prove in a constructive way that multivariate integration in some weighted Sobolev spaces is strongly tractable with -exponent equal to 1, which is the best possible value under a stronger assumption than Sloan and Wozniakowski's assumption. We show that quasi-Monte Carlo algorithms using Niederreiter's -sequences and Sobol sequences achieve the optimal convergence order for any 0$"> independent of the dimension with a worst case deterministic guarantee (where is the number of function evaluations). This implies that strong tractability with the best -exponent can be achieved in appropriate weighted Sobolev spaces by using Niederreiter's -sequences and Sobol sequences.

     

On tractability of weighted integration over bounded and unbounded regions in

We prove that for the space of functions with mixed first derivatives bounded in norm, the weighted integration problem over bounded or unbounded regions is equivalent to the corresponding classical integration problem over the unit cube, provided that the integration domain and weight have product forms. This correspondence yields tractability of the general weighted integration problem.

     

On strong tractability of weighted multivariate integration

We prove that for every dimension and every number of points, there exists a point-set whose -weighted unanchored discrepancy is bounded from above by independently of provided that the sequence has for some (even arbitrarily large) . Here is a positive number that could be chosen arbitrarily close to zero and depends on but not on or . This result yields strong tractability of the corresponding integration problems including approximation of weighted integrals over unbounded domains such as . It also supplements the results that provide an upper bound of the form when .

     

On stability and Bohl exponent of linear singular systems of difference equations with variable coefficients

In this paper we investigate the stability of linear singular systems of difference equations with variable coefficients by the projector-based approach. We study the preservation of uniform/exponential stability when the system coefficients are subject to allowable perturbations. A Bohl–Perron type theorem is obtained which provides a necessary and sufficient condition for the boundedness of solutions of nonhomogenous systems. The notion of Bohl exponent is introduced and we characterize the relation between the exponential stability and the Bohl exponent. Finally, robustness of the Bohl exponent with respect to allowable perturbations is investigated.     

Index-aware model order reduction for differential-algebraic equations

We introduce a model order reduction (MOR) procedure for differential-algebraic equations, which is based on the intrinsic differential equation contained in the starting system and on the remaining algebraic constraints. The decoupling procedure in differential and algebraic part is based on the projector and matrix chain which leads to the definition of tractability index. The differential part can be reduced by using any MOR method, we use Krylov-based projection methods to illustrate our approach. The reduction on the differential part induces a reduction on the algebraic part. In this paper, we present the method for index-1 differential-algebraic equations. We implement numerically this procedure and show numerical evidence of its validity.